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Abstract 

We perform a systematic gradient expansion on kinetic equations and derive the CP- 
violating semiclassical force for fermions propagating in presence of a CP- violating wall at a 
first order electroweak phase transition. The force appears at order h in the flow term of the 
kinetic equation and agrees with the semiclassical force used for baryogenesis computations. 
In particular we consider the force for charginos in both the MSSM and NMSSM. We then 
study the continuity equations for the vector and axial vector currents and stress the role 
of the latter as the one containing baryogenesis sources. We also show that there is no 
CP-violating force for bosons to order % in gradient expansion. 
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1 Introduction 



The creation of a baryon asymmetry at a first-order electroweak phase transition in the 
early universe is an attractive proposal [1] because the elementary particles and interactions 
involved in this process can be tested soon in accelerator experiments. For a successful 
baryogenesis a coalition between CP violation, nonequilibrium thermodynamics and baryon 
number violation is needed [2]. Model calculations require a study of generation and trans- 
port of CP-violating flows arising from interactions of fermions with the expanding phase 
transition fronts. As the problem involves the dynamics of quantum fields in a spatially 
varying background it cannot be treated within the classical transport theory. While fully 
general quantum Boltzmann equations can quite easily be formulated by making use of suit- 
ably truncated Dyson-Schwinger equations for the out-of-equilibrium two-point functions, 
some approximation scheme is needed to derive a set of practically solvable, yet sufficiently 
general transport equations for electroweak baryogenesis (EWBG). 

A fast baryon number violating rate in the unbroken phase is a necessary ingredient of 
any EWBG model. However, to avoid a wash-out of the newly created asymmetry, the 
baryon number violation must turn off in the Higgs phase. As is well known, this is the case 
provided the transition is strong enough [3]. Since for the present experimental bounds on 
the Higgs mass the electroweak phase transition in the Standard Model (SM) is not first- 
order [4] , one is lead to consider extensions of the Standard Model. The natural candidates 
are supersymmetric models, which include the Minimal Supersymmetric Standard Model 
(MSSM) [5] and the nonminimal extension (NMSSM) [6] with an additional Higgs-singlet 
field. These models contain additional scalars which can strengthen the phase transition as 
required for baryogenesis. 

In supersymmetric models the bubble walls are found to be quite slow and thick [7, 8, 9, 
10, 11], in the sense that i> wa n <C c and £ wa n ^> £dB, where £dB ~ ^/T denotes a typical de 
Broglie wave length of thermal particles. The latter condition is of particular importance, 
because it implies that a gradient expansion in terms of £dB/^waii represents a controlled, 
rapidly converging approximation scheme for most of the excitations in the electroweak 
plasma. 
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In the past many heuristic attempts have been made to derive approximate transport 
equations for EWBG [12, 13, 14, 15, 16, 17]. Common to all these methods is the strategy to 
somehow isolate the essential quantum features of the transport in the form of "sources" to be 
inserted into classical transport equations. Baryon production has in this way been computed 
in two doublet models [13, 14, 15, 18], MSSM [16, 19, 20, 21, 22, 23] and NMSSM [24, 25]. 
Different approaches, when applied to the same physical problem, have been found to disagree 
however. In particular the sources from chargino and squark sectors in the MSSM, found 
using WKB-methods [21, 22, 24], are parametrically different from those derived by the use 
of the continuity equations and the relaxation time approximation [19, 20], and by other 
earlier attempts [16, 26]. 

In this paper we present a rigorous first-principle derivation of quantum transport equa- 
tions appropriate for baryogenesis in the limit of thick phase boundaries ignoring collisions. 
We start our analysis by writing the exact Dirac equation of motion for the dynamical Green 
function ( Wightman function) v) = i(Tpp{v)tp a {u)) with a CP- violating spatially vary- 

ing pseudoscalar mass term. For simplicity here we consider particles moving perpendicular 
to the phase boundary, which effectively reduces our problem to 1 + 1 dimensions. The re- 
sults discussed here are not affected in any important way when the general 3+1 dimensional 
case is considered [27]. By performing a Wigner transform we obtain a controlled expansion 
in gradients, or more appropriately, in powers of h. We show that, to first order in H, G < 
admits a spectral decomposition in terms of on-shell quasiparticle excitations. The on-shell 
momenta are set by a dispersion relation derived from the equations of motion and agree 
with the results of [28], where the spectral function A was considered in gradient expan- 
sion. The on-shell distribution functions f s+ and / s _ for particles and antiparticles of spin 
s, respectively, are then shown to obey the following kinetic Liouville equations: 

d t f s ± + v s± d z f s ± + F s ±d k J s ± = 0. (1) 

The quantum information in (1) is entirely contained in the expression for the quasiparticle 
energy ou s ±, which shows up in the expressions for the group velocity v s ± = k z /u s ± and the 
semiclassical force F s ± = u s ±dv s ±/dt, where k z denotes the kinetic momentum. For example, 
we show that in the case of a single chiral fermion moving in a CP-violating background 
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with planar symmetry, represented by a ^-dependent complex mass m(z) 
quasiparticle moving in z-direction with momentum k z has the energy 



m(z)\e ie( - z \ 



a 



uj s ± = ui Q ^h 




(2) 



where u = Jk 2 z + \m\ 2 , and experiences the force 



2u s ± 



m 



s(\m\ 2 0')' 
2l| 



(3) 



We also derive explicit expressions for the semiclassical force for a general case of N mixing 
fermions and in particular for charginos both in the MSSM and NMSSM. We then show 
that in the case of iV mixing bosonic fields, such as the squarks in the (N)MSSM, there is 
no CP-violating force to first order in U. 

Our results agree with recent results obtained by the use of the WKB-approach [22, 24]. 
The WKB-method was originally introduced by Joyce, Prokopec and Turok in [29] and [15] 
and then applied to the MSSM in [21]. The CP-violating velocities and accelerations for 
fermions interacting with a phase transition wall were correctly computed from the WKB- 
dispersion relations in Refs. [29, 15, 21]. The velocity and force in kinetic equations were 
obtained from the Hamilton equations based on canonical momentum. However, when the 
dispersion relation is derived by considering the spectral function in gradient approxima- 
tion [28], the momentum appearing in the Wigner representation is the kinetic momentum. 
The relevance of the kinetic momentum as the true physical variable in the WKB-picture 
was first realized by Cline, Joyce and Kainulainen [22], who also showed how it can be 
consistently incorporated into kinetic theory leading to equations identical to (1-3). 

The outstanding contribution of the present work is a controlled first principle derivation 
of the kinetic equation (1). This is important because of a considerable controversy in 
literature concerning the transport equations to be used for EWBG calculations. Moreover, 
our treatment in principle allows a study of the plasma dynamics beyond first order in H, 
which cannot be achieved by WKB-methods. As an example, at second order in h the full 
equations do not admit the spectral decomposition solution for G < ; for a scalar field this 
problem is considered in [30]. 
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Let us mention that in a related work [31] the Liouville equations for fermions in presence 
of a classical gauge field have been considered in gradient approximation. The crucial role of 
the constraint equations in the derivation of the kinetic equations was then stressed in [32, 
33]. The problem of a pseudo-scalar mass term in kinetic equations has been considered in 
Refs. [34, 32], but these authors discussed the flow term only to zeroth (classical) order in 
h, whereas the spin dependent force essential for EWBG arises only at quantum level, as we 
show here. 

Inclusion of interaction terms gives rise to yet another source which is of first order in 
h, namely the spontaneous baryogenesis (SBG) source of Ref. [35]. The SBG source appears 
because the CP-violating split in the dispersion relation causes the CP-conjugate states to 
relax towards different local equilibria in the bubble wall. Thus a first principle derivation 
of the SBG source requires not only a consistent expansion in h, which we have done, but 
also a consistent expansion in relevant coupling constants. The latter is necessarily a model 
dependent problem and shall be considered elsewhere. However, to facilitate comparison 
with literature we derive the vector current to order H in gradient expansion which displays 
the SBG source in the relaxation time approximation. When applied to the MSSM, our 
results differ from Refs. [36], [17] and [20]. 

The paper is organized as follows. In section 2 we derive the Liouville equations for a 
single Dirac fermion with a spatially varying complex mass term. In section 3 we generalize 
these results to the case of N mixing fermionic fields and study the case of mixing charginos 
in both the MSSM and NMSSM. We then in section 4 consider the case of N mixing scalar 
fields and show that there is no CP-violating source to first order in H. In section 5 we 
study the continuity equations for both vector and axial vector current, and spontaneous 
baryogenesis in the relaxation time approximation, and make a comparison with literature. 
For example, in contrast to what is claimed in [37] and [17], we find that the continuity 
equation for the vector current contains no CP-violating source in the absence of collisions. 
On the other hand, in the continuity equation for the axial current there are CP-violating 
sources that can be related to higher moments expansion of the semiclassical Boltzmann 
equation. Finally, section 6 contains a discussion and summary. 
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2 Fermionic field with a complex mass 



We first consider the dynamics of a fermionic field with a complex spatially varying mass 
term. More precisely, we take our system to be described by the effective lagrangian 

C = ii)pi) - 4> L mip R - i> R m*^ L + C mt , (4) 

where £ int contains interactions and 

m(x) = m R (x) + im I (x) = \m(x)\e ie ^ (5) 

is a space-time dependent mass term arising from an interaction with some CP-violating 
scalar field condensate. We are primarily interested in the case where m arises from the 
Higgs field condensate at a first order electroweak phase transition. As the bubbles of broken 
phase grow several orders of magnitude larger than the wall width before coalescence, the 
wall can be approximated by a planar interface to good accuracy. We therefore consider 
a mass term in the bubble wall frame which is only a function of the spatial coordinate 
orthogonal to the wall, m = m(z). 

Our focus in this paper is on the semiclassical treatment of fermions in presence of 
background fields. In particular we derive the flow term for a fermionic kinetic equation with 
a nontrivial force induced by the CP-violating mass parameter (5). To keep the discussion 
simple, we relegate the explicit treatment of collision terms and self-energy corrections which 
are induced by the specific interactions included in £i nt , to later publications [27]. Moreover, 
we are interested in the case of wide walls, so that the de Broglie wave length £^b of a typical 
excitation is small in comparison with the wall width, £^b £ W aii- This condition is amply 
satisfied at the electroweak phase transition, where typically £<ib ~ 1/T and £ wa u ~ 10/T [8]. 

With the above assumptions we now develop the equations of motion for the Wightman 
function 

Gap(u, v) = % (^(v)^(m)) (6) 

in a consistent expansion in gradients of the background fields. Here (•) denotes the ex- 
pectation value with respect to the initial state. The function G K describes the statistical 
properties of an out-of-equilibrium system. It corresponds to the off-diagonal part of the 
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fermionic two-point function in the Schwinger-Keldysh formalism [38], which indeed is the 
method of choice to derive the equations of motion for G < including interactions. However, 
in their absence all one needs is the familiar Dirac equation. Dropping interactions we have 
from Eq. (4): 

(i $ u -m R (u) - rfm^u)^ ip(u) = . (7) 
Multiplying (7) from the left by the spinor iip(v) and taking the expectation value one finds: 

(ifi u -m R (u)-i 1 5 m I (u))G < (u,v) =0. (8) 

This equation and the hermiticity property 



t 



ii s G < {u,v)\ = i^G < (v,u), (9) 

which can be immediately inferred from the definition (6), completely specifies G < . 

In order to study equation (8) in gradient expansion we perform a Wigner transform 
of G < to the mixed representation, i.e. a Fourier transform with respect to the relative 
coordinate r = u — v: 

G < (x,k) = f d 4 re ikr G < (x + r/2,x-r/2), (10) 

where x = (u + v)/2 denotes the center-of-mass coordinate. The crucial advantage of the 
representation (10) is that it separates the internal fluctuation scales, described by momenta 
k, from the external ones which show up as a dependence of G < on x, and thus gives us the 
chance of exploiting possible hierarchies between these scales. In the Wigner representation 
equation (8) becomes 

(jfc - m {x) - im 5 (x)-f 5 ) G< = 0, (11) 
where we use the following convenient shorthand notation 

kfj, = k^+'-d^ (12) 
7710(5) = m R(I) e-^- dk . (13) 

The original local equation (8) for G < is thus transformed into an equation involving an 
infinite series in gradients. This in fact can be viewed as an expansion in powers of the 
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Planck constant h. We have set h — > 1, but a dimensionful h can at any stage be easily 
restored by the simple replacements d x ^ hd x and G K — > /i -1 G < . 

Because of the planar symmetry (here we do not consider initial states of the plasma 
that break this symmetry) G < can depend only on the spatial coordinate orthogonal to the 
wall, z = x 3 . We also consider equation (11) only in a frame where the momentum parallel 
to the wall vanishes, k\\ = 0. This last assumption effectively reduces the problem to 1+1 
dimensions, and we can cast equation (11) into the form 

(k + kl°l 3 ~ rh ol ° + im 57 7 5 ) i 1 Q G < = 0, (14) 

where k = k + ^d t and k z = k z — ^d z . The differential operator in (14) is entirely spanned by 
a closed 1+1-dimensional subalgebra of the full 3+1-dimensional Clifford algebra. Moreover, 
it commutes with the operator S 3 = 7°7 3 7 5 , which measures spin s in z-direction. s is thus 
a good quantum number in the frame k\\ — 0, which motivates to seek solutions for i / y°G < 
which satisfy S 3 i r y°Gf = i^GfS 3 = sij°Gf. Working in a convenient chiral representation 
this condition leads immediately to the following spinor structure: 

- tl °Gf = ±(l + sa 3 )®gf, (15) 

where a 3 is the Pauli matrix referring to spin in z-direction. gf has indices in the remaining 
two dimensional chiral space and can also be written in terms of the Pauli matrices p l as 
follows: 

9f = \{^ + ^j)- (16) 
The decomposition (15) contains the implicit assumption that i^ Q G < does not mix spins. 
We expect that this is a good approximation in the electroweak plasma. The signs and 
normalizations in (15-16) are chosen such that g$ measures the number density of particles 
with spin s in phase space. With these simplifications we can now reduce our original 4x4 
problem to a two-dimensional one by effecting the replacements 

7° - P\ -i7°7 5 - P 2 , -7 5 - P 3 (17) 

and we find 

(k Q - sk z p 3 - p 1 m - p 2 m^j gf = 0. (18) 



A similar procedure is used in [28] for a treatment of the fermionic propagator. Equation (18) 
may look simple, but it still constitutes a set of four coupled complex (or eight coupled 
real) differential equations, which we shall now analyze. We first find the four independent 
complex equations for g s a (a = 0,i) by multiplying (18) successively by 1 and p\ and taking 
the trace: 

k g'o - sk z g s 3 - m Q g{ - m 5 g% = (19) 
ko9t - sk z g s Q - im g s 2 + irh b g{ = (20) 
k Q g{ + isk z g s 2 - m g-Q - im 5 g s 3 = (21) 
k g s 2 - isk z g{ + im g s 3 - m^g s Q = 0. (22) 

As a consequence of Eq. (9) the matrices gf are hermitean so that g* are real functions. 
We then have twice as many equations as independent functions corresponding to real and 
imaginary parts of Eqs. (19-22), and hence one half of the equations must correspond to the 
constraints on the solutions of the other half; those equations are kinetic equations. This 
was first pointed out in the context of kinetics of fermions by Zhuang and Heinz [32]. As one 
sees from (12), no time derivatives appear in the real parts of (19-22) so that they indeed 
provide four constraint equations (CE) on the solutions of four kinetic equations (KE). These 
contain time derivatives and correspond to the imaginary parts of (19-22). 

Because we have put no restrictions to the form of the m-operators, Eqs. (19-22) are still 

— * 

valid to any order in gradients in the frame where k\\ =0. In what follows we assume that 
the mass is a slowly varying function of x and truncate gradient expansion at second order, 
which is the lowest order at which CP-violating effects can be discussed consistently. This 
method is not adequate for problems involving quantum mechanical reflection, which require 
nonperturbative treatment in H. With the truncation, m (5) in (19-22) simplifies to 

i 1 

m ( 5 ) ^ m m + -m' R(I) d kz - -m" R{I) d 2 kz . (23) 

Even with this truncation, we are facing a problem involving eight coupled second order 
partial differential equations. Our next task is to reduce these to a single equation governing 
the dynamics of the fermionic two-point function. 
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2.1 Constraint equations 

Let us first consider the constraint equations. They consist of four homogeneous equations 
for four functions, which implies that there is one constraint which gives rise to the dispersion 
relation. While this property remains true to any order in gradients (or equivalently in H), 
we only need to work to first order to find the nontrivial result we are looking for. To this 
order we have 

ko9o ~ sk z9 s 3 ~ m R9i ~ m i92 = (24) 
ko93 ~ sk z9o + 2 m> R d kz92 ~ l^n'idk z g[ = (25) 

Mi + \dz92 ~ m R g s + \vn!jd kz gl = (26) 

s 1 
ko9 s 2 ~ ^ d z9 S i - m ig s - -m' R d kz g s 3 = 0. (27) 

We first use the constraint equations (25-27) iteratively up to first order to express gf, g s 2 
and g\ in terms of g^ and d kz g^ and then insert the results into (24). Remarkably all 
terms proportional to <9fc z (?o cancel and we find that to first order in gradients g$ satisfies the 
algebraic equation 

^-A; 2 2 -|m| 2 + ^|m| 2 ^)^ = 0. (28) 
This admits the spectral solution g$ = 7rn s \k \5(Ql), which can also be written as 

9o = Y,7rj—ns6(koTu s ±). (29) 

Here n s (k , k z , z) are nonsingular functions which are, as we show below, related to the on- 
shell distribution functions. The indices ± refer to the sign of k , to be eventually related to 
particles and antiparticles. The energy u s ± is specified by the roots of the equation 

Q 2 S = kl-k 2 z - |m| 2 + f \m\ 2 6' = 0, (30) 

and the normalization factor Z s ± is defined as 

Z s ± = ~ \dk ti 2 s \k =±oj s ±- (31) 

Zuj s ± 
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To first order in gradients these can be solved iteratively: 



Z s ± = Its 



u s ± = too =F s 




u = 



ml 2 



(32) 



(33) 



Equation (32) defines the physical dispersion relations for particles and antiparticles of a 
given spin s. Due to the derivative corrections the spin degeneracy is lifted at first order in 
gradients and hence particles (antiparticles) of different spin experience different accelera- 
tions in a spatially varying background, as we shall see in more detail below. 

Solution (29) nicely illustrates how the constraint equations operate. Solving (24-27) 
consistently to first order accuracy constrains the solutions of the kinetic equation to sharp, 
locally varying energy shells given by (32). One way of understanding this is as follows. 
The CP-violating phase 6' in (28) can be related to an axial gauge field [29, 15] which, to 
leading order in gradients, lifts the degeneracy in the dispersion relation, but does not spoil 
the quasiparticle picture, just as it is the case with a vector gauge field. We should note 
however, that the confinement to sharp energy shells does not persist beyond first order in 
gradients. While for noninteracting fermions one can always express <7i 2 3 i n terms of g$, 
at higher orders more complicated derivative structures arise, as the constraint cannot be 
written as a simple algebraic equation with a spectral solution. For a treatment of such a 
situation in the case of a scalar field see [30]. 

Eq. (32) is identical with the results derived earlier by WKB-methods in [22] and simul- 
taneously via the field-theoretic technique of spectral integrals in [28] . From the WKB-point 
of view the present derivation comes as a welcome verification of the result obtained by an 
intuitive, but less fundamental approach. The agreement with the field theoretical calcula- 
tion of [28] on the other hand was to be expected. In [28] it was shown that the integral 
over the spectral function, defined as a difference of the retarded and advanced Green func- 
tions A = (i/2)(G ret — G a dv), projects test functions onto energy shells (32); however, in 
the collisionless limit A satisfies the same equation of motion (8) as the Wightman function 
G < . The trace of rfA in particular satisfies equation (28), and can be obtained from (29) 
by the replacement n s — > 1. We can hence immediately check the sum-rule to the accuracy 
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at which we are working. Indeed, 

r —Tr^As = £ f ^^o^ s± ) = 1. (34) 

We finally note that equation (30) has additional poles at ko ~ s|m| 2 #'/2u;o, which we 
have left out in the decomposition (29). These poles correspond to unphysical, but harmless, 
tachyonic modes which arise only because our solutions for the constraint equations involve 
an expansion in inverse powers of ko, which breaks down already for ko much larger than the 
value associated with these poles. Note that their contribution to the sum rule (34) vanishes 
when summed over spins. 

2.2 Kinetic equations 

We now turn our attention to the kinetic equations. We are primarily interested in the 
equation for which carries information on the particle density in phase space. From (19) 
we have 

d t 9o + sd z9 s 3 - m 'ndk z g s i ~ m'jd kz g s 2 = 0, (35) 

which is correct up to second order in gradients (first order in H). Just as in the previous 
section we use the constraint equations (25-27) to express g(, g s 2 and g\ in terms of g^ and 
arrive at an equation for g$ alone. To second order in gradients (first order in K) it reads 

k d t g s + k z d z g s - fym\ 2 ' - ^-(|m| V)') d kz g s = 0. (36) 

We have so far used three out of four constraint equations. To find the acceptable solu- 
tions satisfying all constraints, we must yet impose the restriction onto the functional space 
spanned by (29). Because of the 5-function in the decomposition (29) this is of course trivial. 
All we need to do is to insert (29) into (36) and integrate over the positive and negative 
frequencies k . We then get the following form for the Liouville equation 

d t f s ± + v s ±d z f s ± + F s ±d k J s ± = 0, (37) 

where 

fs+ = n s (uj s+ , k z , z) 

/«_ = l-n 8 (-w a -,-k x ,z) (38) 
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are the distribution functions for particles and antiparticles with spin s, respectively. These 
definitions are motivated by the equilibrium result n e s q = l/(e /3fc ° + 1), where (3 = 1/T is the 
inverse temperature. The quasiparticle group velocity v s ± appearing in Eq. (37) is given by 

v s ± = (39) 

where k z is the kinetic momentum. The spin-dependent and CP-violating semiclassical force 
reads 

f _ ±= _hr ± ^ (40) 

2uj s ± 2cjq 

Equations (37-40) are among the main results of this paper. Incidentally, the form (40) 
for the semiclassical force F s ± = u) s ±dv s ±/ dt was already found by Joyce, Prokopec and 
Turok [29]. To obtain kinetic equations by WKB-methods the authors of [15] used canonical 
variables however. The resulting equations are not invariant under reparametrization of 
the wave functions, and hence care is required when specifying local thermal equilibrium in 
derivation of transport equations relevant for baryogenesis. Cline, Joyce and Kainulainen [22] 
introduced the kinetic momentum as a physical variable in the kinetic equations and obtained 
the unique reparametrization invariant transport equations identical with (37) and (39-40). 
The outstanding contribution of the present work is in a controlled first principle derivation 
of these equations without any a priori assumptions. 

Let us finally note that equation (36) could have been obtained by taking the bilinear 
O-derivative of the constraint equation (28), where the O-derivative is defined by 

0{a}{b} = i (d t a d ko b - d z a d kz b - d ko a d t b + d kz a d z b) . (41) 

This is no coincidence, and even more generally, in the collisionless limit the kinetic equation 
can be obtained by effecting the tan O-derivative on the constraint equation. 

2.3 Currents 

It is instructive to study the expressions for physical currents in order to shed light on various 
functions we have encountered in our derivation. Of particular relevance for baryogenesis 
are the vector and axial vector currents. By making use of (6) and (15-16) one finds 

d 2 k 



s=±l ■ 
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# = (^tViK*)) = Zj — i&sgl), (42) 

where we have restricted ourselves to 1+1- dimensions so that d 2 k = dk z dko- This shows 
that Qq is the usual number density in phase space, whereas gj represents the axial charge 
density. An important consequence of the constraint equations is that there is only one 
independent dynamical function, here chosen to be while all others can be related to g$ 
via the constraint equations (25-27). In particular, (?| can be written as 

The nontrivial gradient correction appears as a total derivative and hence vanishes upon the 
/^-integration. Using the decomposition (29) one finds 

& = l£k^'*= (44) 

oL± = s f (v s ±, l) / s± = s s k* I (!> 7—) /«± > ( 45 ) 

where denotes the sign of fc z , and we discarded the vacuum contribution. In the last 
step we used dk z 0J s ± = Z~±k z jio s ±. The lower limit in the tu-integrals is \m\ =F s9'/2. The 
functions f s ± are the correctly normalized distribution functions, and they retain the correct 
physical interpretation in that the particle flux is not affected by CP-violating effects, while 
the density may be either enhanced or suppressed, as given by the inverse velocity. The 
current (44) was computed by WKB-methods in [22]. The correct result for j®± was also 
found in [28] by field-theoretical methods. However, j%± found in [28] does not agree with 
(44) because the spinor structure used for i^y°G < was too simple. 

2.4 Interactions 

In all discussions above, we have left out the effects of interactions. This was done to avoid 
the necessity to use the full machinery of the Schwinger-Keldysh formalism [38], and to 
keep things as simple as possible. Moreover, unlike the treatment of the flow term of the 
kinetic equation presented above, including interactions is necessarily a model dependent 
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task. Nevertheless it is a simple matter to write a formally exact equation of motion for G < 
including the collision terms. Instead of (8) we then have 

(ip-m R -ima 5 )iG < -S R 0?G < = £< G R + i(s> G< - G< £>), (46) 

where A B{u,v) = / g?k; A(u, w)B(w, v). Gr is the real part of the (retarded) propa- 
gator, and the function £r contains the real part of the self-energy corrections including 
the singular (tadpole) interactions which can be resummed to a renormalized mass term. 
The terms in parentheses give rise to the usual collision term, where the self-energies 
arise only from nonlocal loop contributions to the Dyson-Schwinger equations. However, as 
mentioned above, the exact form of these terms depends on the theory considered. While 
their treatment is not conceptually difficult, their inclusion brings a considerable amount of 
technical complications. We shall consider the problem of including collisions elsewhere [27]. 

3 Mixing fermionic fields 

In practical applications, such as in supersymmetric models, one needs to consider cases 
where several fermion flavours are mixed by a spatially varying mass matrix. We therefore 
consider a theory with the mass lagrangian 

Anass = -^LM^R - $ R M^ L , (47) 

where M is a complex (in general nonhermitean) N x N matrix with spatially varying 
components. We denote the flavour degree of freedom by an additional index i to the spinor 
ip a>i (x), so the Wightman function becomes a matrix in the product space of spinor and 
flavour: 

Ga/3,ij( U , V ) = * (4>/3,j(v)ll>a,i( U )) ■ ( 48 ) 

The flavour degree of freedom plays no role in the derivation of the equations of motion for 
Gap,ij i n t ne steps analogous to going from (8) to (19-22) in the single fermion field case, 
because those steps dealt only with the spinor structure of G < . We can thus immediately 
write an equation analogous to (18): 

(k - sk z p 3 - p l Mo - P 2 M 5 ) gf = 0. (49) 
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The sole, but significant difference to (18) is that gf is now an N x iV-matrix in the flavour 
space, and the mass terms have become N x iV-matrix operators 

M = \{M + M ] ) (50) 
M 5 = -l(M-Mt) (51) 

where M = Me^ dz&k ^ and = M^ei dz9k * . The extra flavour structure of course complicates 
the solution of equation (49) and it turns out to be convenient to perform a rotation to 
the basis where the lowest order mass matrix is diagonal. Because M in general can be 
nonhermitean, the diagonalization requires a biunitary transformation 

UMV^ = M d , (52) 

where U and V are the unitary matrices which diagonalize the hermitean matrices MM^ 
and M'M, respectively. After the rotation we can write (49) in the component form in the 
diagonal basis as 

(k + l -V-)g s od - s ( kz - l -V;)g s 3d - M od g s ld - M bd g s 2d = (53) 
{ko+ l -V;)g s M - s(k z - l -V;)g s 0d -iM 0d g s 2d + iM 5d g s ld = (54) 
(ho + % -Vt)g{ d + is(k z - % -Vt)g s 2d - M Gd g s 0d - iM 5d g s 3d = (55) 
(ko + ^t)g S 2 d - is(k z - \^ > T)9\d + i] ^od93d ~ M 5d g s 0d = 0. (56) 
The 'covariant derivatives' appearing in (53-56) are defined as 

V f = dt-i[E t , -}--is[A z , -] ± (57) 
Vf = d z -i[Z z ,-}_-is[A t ,-} ± (58) 

where the brackets [•,•]- refer to commutators and [•,•]+ to anticommutators and 

S M = l -(yd^ + Ud^) (59) 
A M = l -(Vd^ -Ud^). (60) 
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It should be noted that, while the relation (52) allows arbitrary phase redefinitions U — > wU 
and V — > wV, where w is any diagonal matrix with \wa\ = 1, the operator remains 
invariant under these transformations. This reparametrization freedom is exactly what leads 
to the apparent 'gauge' dependence of the results in the WKB-approach [22]. As we show 
below, only the diagonal elements of A M contribute to the constraint and kinetic equations 
to order H, which then implies that our results are reparametrization invariant. Finally, the 
new mass operators in the diagonal basis are given by 

M Qd = ^(UMV^ + VMW) (61) 
M bd = --(UMV ] -VM ] U ] ) . (62) 

The constraint and kinetic equations now correspond to the hermitean and antihermitean 
parts of (53-56), respectively. Because of the matrix structure, these equations contain a 
number of commutator and anticommutator terms involving g* and the various matrix- 
operators. However, we shall now argue that (53-56) can effectively be taken to be diagonal 
to the order at which we are working. Indeed, in the propagating basis (52) the off-diagonal 
terms are obviously suppressed by h when compared to the diagonal elements. On the 
other hand, they appear in the diagonal equations through the commutators h~ 1 [hT, z , g^] = 
[T, z ,g ad ] and [A z ,g ad ], and thus at the same order as the diagonal elements. This then 
immediately implies that, when the dynamics of CP- violating densities is considered, the off- 
diagonals contribute at second order in h in the diagonal equations and can be consistently 
neglected. With this it is now straightforward to show that, to first order accuracy in h, the 
constraint equations reduce to the following equations for the diagonal entries of g s ad : 

k A - sk z9 s 3 d ~ m R9ld ~ m i92d = ( 63 ) 
kog s 3d - sk z g s od + -m' R d kz g s 2d - -m'jd^g^ = (64) 

s 1 

kogla + 2 9 *0m + sA zd9t d - m R g s od + -fh'jd kz g s zd = (65) 

s 1 

k o9 S 2d - 2 dz 9td + sAzdg S 2d - m i9od - 2"4 9 *.0M = °> ( 66 ) 

where m^i are the real and imaginary parts of the eigenvalues of Md, A z d is the diagonal 
part of A z and fh' R = m' R — 2miA z d and fh'j = w! l + 1m^A z d. Similarly, the kinetic equation 



16 



for t° first order in h becomes 

dtg s od + sd z g s 3d - m' R d kz g[ d - mjd kz g & 2d = 0. (67) 

Following our treatment in section 2, it is now straightforward to eliminate g s ad (a = 1,2,3) 
from equations (63) and (67) to obtain the constraint equation 

(k 2 - k\ - |M d | 2 + ^|M,| 2 0') g° d = (68) 

and the kinetic equation 

k d t g s 0d + k z d z g s 0d - (\\M d \ 2 ' - ^-{\M d \ 2 Q')')d kz g^ d = (69) 

for the number density function g$ d alone. Eqs. (68) and (69) are analogous to the one 
field equations (28) and (36). The difference is that here we have N different equations 
corresponding to N diagonal elements of g$ d in the mass eigenbasis. The derivative of the 
effective angle 9' appearing in (68-69) is defined as 

Q' = Q d + 2A zd (70) 

where the angles are the complex phases of the elements in the diagonal mass matrix: 
M d = \M d \e i0d . Since the mixing field contribution shows up as a shift in the derivative of the 
pseudoscalar phase 2A zd , it implies that equations (68-69) are manifestly reparametrization 
invariant. We can convert (70) to an alternative form in terms of the original mass matrix 
and the rotation matrix U : 

|M d | 2 9' = -ilm(f/(MM /t - M'M^)U ] ) d , (71) 

which is also manifestly reparametrization invariant. This expression will be convenient for 
discussion of the chargino sector of the MSSM and the NMSSM in sections 3.1 and 3.2. 

The final steps in going from equations (68) and (69) to kinetic equations for the mass 
eigenmodes are exactly analogous to the single fermionic field case: each mass eigenmode 
i gets projected to its own energy shell ou S i± given by (68), and the corresponding spectral 
decomposition density function f si ± obeys a semiclassical Boltzmann equation identical to 
(37) 

dtf si ± + v si± d z f si± + F si ±d k J si± = , (72) 
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with the corresponding group velocity 



v.si± = (73) 



Fst± = -L-^± Vl * 2 (74) 



and the CP-violating semiclassical force 

j si ± 2uof 

computed from expression (70) or (71). It has been shown in [29, 15, 22], that the spin- 
dependent term in F si ± gives rise to a CP-violating source proportional to |Mj| 2 6^ in the 
diffusion equations. We can therefore loosely call this factor the 'source', and proceed to 
compute it in some special cases. 

3.1 Charginos in the MSSM 

We first compute the source in the transport equations for charginos in the MSSM. The 
chargino mass term reads 

* fl M* L + h.c, (75) 

where = (W^,hf R ) T and ^/l = ^Wli^2,l) T are the chiral fields in the basis of winos. 
The mass matrix reads 

where H 1 and H 2 are the Higgs field vacuum expectation values and /i and m 2 are the 
soft supersymmetry breaking parameters. For a realistic choice of parameters there is no 
spontaneous CP-violation in the MSSM, so to a good approximation we can take the Higgs 
vev's to be real [25, 23]. The matrix U in (52) can be parametrized as [22] 

^ H(A + A) a \ 



with 



a = g{m 2 H 1 + n*Hl) (78) 
A = \ m ^-\^ + g\hl-h\) (79) 
A = v /A 2 + 4|a| 2 , (80) 
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where hi = \Hi\ are normalized such that the tree level W^-boson mass is = g 2 h 2 /2 : 
h 2 = h\ + h\. The physical chargino mass eigenvalues are given by 

ml = \(\m 2 \ 2 +\tf + g 2 h 2 )±±. (81) 

Upon inserting (76) and (77) into (71) it is straightforward to show that the source term for 
charginos becomes 

rnle' ± = ^Z(Lim 2 )(h 1 h 2 y, (82) 

where + (0_) corresponds to the higgsino-like state when > |m 2 | < |m 2 |). CP 
violation is here mediated via the parameters /i, m 2 and may in fact be large [39]. The 
result (82) is in perfect agreement with the chargino source obtained by a WKB method in 
Ref. [22]. 

3.2 Charginos in the NMSSM 

In the NMSSM there is an additional singlet field S in the Higgs sector. The singlet field 
couples to higgsinos, and hence the higgsino-higgsino component in the chargino mass matrix 
(76) is generalized in the NMSSM: 

fi — > fi = fi + XS, (83) 

where A is the coupling for higgs(ino)-higgs(ino)-singlet interaction. Another consequence 
of this extension is the possibility to have spontaneous transitional CP-violation [25], so the 
Higgs fields Hi are in general complex. When the parameters a, A and A are defined as in 
equation (80) with /i — > jl, the matrix U in equation (77) still diagonalizes MM^ . 

In the NMSSM we must account for the complex dynamical phases of the higgs fields, 
and hence we have to be more careful with our definition of the higgs doublets. Our choice 
of writing the mass matrices (76) corresponds to parametrizing the higgs doublets $j as [40]: 

•>-(*£)■ M 

where hf are the charged higgs fields. Only one of the higgs phases 9i is physical, while the 
other gets eaten by the gauge fields in the unitary gauge. We wish to choose the physical 
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phase in such a way that the corresponding field does not couple to the neutral weak boson. 
Given the parametrization (84), this condition implies that 



h{9[ = hie!,. (85) 
Using the gauge constraint (85) we can write 

e[ = f 2 e>, 9> = §9>, (86) 

where h 2 = h\ + h 2 , and 9 = Q\ + 9 2 is the physical CP- violating phase. 

To get the explicit form for the CP-violating term, we insert the NMSSM mass matrix (76) 
into (71). After some algebra one finds the following three terms giving rise to CP-violating 
sources in the NMSSM: 

©NMSSM = ©fcife + + ©s > ( 87 ) 
The first term is the following generalization of the chargino source (82): 

™I©U± = ( W ' ( 88 ) 



A 

for the case involving a new scalar field S and possibly complex higgs fields. However, there 
are two new types of terms in the NMSSM. The term Q' e is proportional to a derivative of 
the CP-violating phase 9 in the Higgs sector, and reads 

m%& e± = ( (A ± (|m 2 | 2 + m)^T mim^hhi) . (89) 



Finally, the source Q' s can be written as a derivative of the singlet condensate: 

<®'s± = ±^Q(m 2 H 1 H 2 S') + ^(A±m 2 + g 2 h 2 - \m 2 \ 2 ))^*S') . (90) 

In all formulae (88-90) the mass eigenvalues m\ can be read off from equation (81) with the 
replacement /i — > jl, where jl is given by Eq. (83). Baryogenesis in the NMSSM from the 
semiclassical force has been studied in Ref. [24]. 
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4 Mixing bosonic fields 



Here we first show that, unlike for fermions, the constraint and kinetic equations for mixing 
bosons acquire no gradient correction to first order in h in the collisionless limit. We then 
derive the constraint and kinetic equations accurate to second order in gradients which can 
be used as a starting point for baryogenesis calculations. N mixing bosonic fields with a 
spatially varying mass matrix M 2 obey the Klein-Gordon equation 

(□ U + M 2 H)0( M )=O, (91) 

where is an N- dimensional vector whose components are coupled by the hermitean mass 
matrix M 2 . Multiplying (91) from the left by —i(ft(y) and taking the expectation value with 
respect to the initial state we get 

(n u + M 2 {u))G < {u,v) = 0. (92) 
where the Wightman function G < is defined as 

G < (u,v) = -i(<i>\v)<i>(u)). (93) 
After performing the Wigner transform, equation (92) becomes 

Qd 2 _ k 2 - ik ■ d + M 2 e~5 9 ■ 9k^ G < = Q ( 94 ) 

In the case when TV = 1 it is immediately clear that the first quantum correction to the 
constraint equation (the real part of (94)) is of second order and to the kinetic equation 
(imaginary part) of third order in gradients (second order in h). To extract the spectral 
information to second order in Ti is quite delicate since the constraint equation in (94) 
contains derivatives [30]. In the case of more than one mixing fields it is convenient to 
rotate into the mass eigenbasis, just as in the fermionic case: 

M 2 = UM 2 U\ (95) 

where U is a unitary matrix. In the propagating basis the equation (94) becomes 

(bp 2 -k 2 -ik-V + M 2 e-^- 9 ^ G< = 0, (96) 
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where G^ = UG < W and the 'covariant' derivative is defined as: 

P„ = ^-i[S M , •], H M = iUd^U^. (97) 

Since (G^)^ = — and 2^ = £> M , we identify the antihermitean part of (96) as the con- 
straint equation: 

- 2k 2 G< + [Ml + -V\ Gf}-i[k-V + Ml G<] = 0, (98) 
and the hermitean part is the kinetic equation 

{k.V + MlG<}- 

where we defined 



M c 2 + X -V\G< 



0, (99) 



Ml = Ml cos^V-d k 



Ml = Mlsm^V-d k . (100) 

We now use the analogous argument as in the fermionic case in section 3. The off-diagonal 
elements of G^ in (98-99) are sourced by the diagonal elements through the terms involving 
commutators which are suppressed by at least H with respect to the diagonal elements. This 
implies that, in order to capture the leading order nontrivial effect in gradients, we can 
work in the diagonal (semiclassical) approximation for G^. By inspection of (98-99) we can 
now immediately write the constraint and kinetic equations in the diagonal approximation 
accurate to order H as follows 

(k 2 -Ml)G< = (101) 
(k-d+ l -(dMl)- dk) G< = 0. (102) 

In contrast to the fermionic equivalent (68), these equations contain no CP- violating correc- 
tions to order h and display only the usual classical CP-conserving term associated with the 
mass eigenvalues. This analysis is relevant for example for calculation of the CP-violating 
force in the stop sector q = (t^,t^) T of the MSSM, in which the mass matrix reads 

M 2_( m% y{A*H 2 + ^H l )\ 
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where ttiq and denote the sum of the soft susy-breaking masses, including D-terms and 
m 2 = y 2 H\. Our analysis immediately implies that for squarks there is no CP-violating 
correction to the dispersion relation at first order in gradients, and hence there is no CP- 
violating semiclassical force in the kinetic equation at order H. 

5 Continuity equations and CP- violating sources 

The quantities eventually relevant for baryogenesis are CP-violating fluxes. We now investi- 
gate how the CP-violating sources appear in the equations for the divergences of the vector 
and axial vector currents. Let us first consider the divergence of the vector current. We have 

d,f = d^x^ix)) (104) 

= /^W-^V) (105) 

where the derivative in the integral expression is taken with respect to the center-of-mass 
coordinate. Using the decomposition (16) and the constraint equations to write g^ in terms 
of g^ (i is the flavour index), Eq. (16) is easily shown to reduce to just a momentum integral 
over the kinetic equation (69). Hence we get the usual continuity equation 

d/J&± = d t n si ± + d z (n si ±u si ±) = 0, (106) 

showing that the vector current is conserved, and contains no sources. The fluid density n S i± 
and velocity u si ± = (v si ±) are defined as 



r d 2 k 

i: = y + (2vp^ 

d 2 k ( k v " 



«*«•> - (107) 

where J + = J k >0 denotes integration over the positive frequencies. The fluid density n si ± 
should not be confused with the phase space density n s in (38). To get the density and 
velocity moments for antiparticles one should integrate over the negative frequencies and 
make use of (38). Here we have again restricted ourselves to 1+1 dimensions; in 3+1 di- 
mensions the expressions differ in detail, but not in essence [27]. No source appears in (106) 
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simply because the semiclassical force term in (69) reduces to vanishing boundary terms at 
k z — > ±00. In the case of a single Dirac fermion this result remains valid to any order in 
gradients. In a general case of mixing fields the dynamics of off-diagonal elements of g$ may 
induce sources, which would however be higher than first order in h. Our proof that there are 
no CP-violating sources to the continuity equation for the vector current is contrary to the 
results of Refs. [37, 17]. Our derivation is more general than that of [37] in that by treating 
mass as part of the flow term it includes the "mass resummation" of Refs. [37, 17] to infinite 
order, it is not based on any particular Ansatz for g^ d , and finally, but most importantly, we 
took correct account of the constraint equations. The fact that we have not treated collisions 
terms here does not resolve the differences, because the collisional contributions arise only 
from nonsingular loop diagrams which were not treated in Ref. [37] either. 

We next consider the continuity equation for the axial vector current. This can be 
obtained from Eq. (106) by replacing 7 M by 7 M 7 5 . The presence of 7 s essentially changes 
the roles of g% and g$, so that the axial divergence reduces to an integral over the kinetic 
equation for g|, which can be inferred from (54): 

dtgl + sd z g s + 2 [Mj - \M'id 2 k }j g{ - 2 [m r - \m h r ^ g s 2 = 0. (108) 

It is thus easy to see that the divergence of the axial current acquires the expected form: 

d^ s = -2iM R $ sl 5 ^ s ) - 2M J (^ S >, (109) 

where (4> s ips) = Ik 9i denotes the scalar, and {4> s , y 5 ip s ) — i J k g 2 the pseudoscalar density, and 
J k = J d 2 k/(2-K) 2 . We now make use of the constraint equations (53-56) to express g[, g 2 and 
(?| in (108) in terms of g s Q (to second order in gradients), diagonalize (108) by rotating into 
the propagating basis, where we can take g^ d to be diagonal to order H accuracy, and finally 
integrate over the momenta. We thus arrive at the following equation for the axial current 
divergence: 

d t (n si ±u si ±) + d z n si ± - [\Mi\ 2 d z + ^\M i \ 2 ^jX 2s i± 

± s (jm,| 2 ©^ + ^flMife;)') x 3si± = 0, (no) 
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where we defined 

As usual, one has to introduce some truncation scheme to close the equations to two unknown 
quantities (here n si ± and u S i±). There is of course some freedom as to how to do this step, 
and one might implement the truncation in (110) by replacing the distributions / S j± in 
the X ps i±-integrals by the equilibrium ones. However, it is instructive to observe that, by 
extracting a total derivative from the |Mj| 2 -terms in (110), the corresponding integrals can 
be combined with the <9 z n si ±-term to give a second velocity moment term, in terms of which 
(110) becomes simply 

d t (n si ±u si ±) + d g (n si ±(v 2 si± )^ = S si ±, (112) 
where the source S si ± is given by the average over the semiclassical force (74) divided by k : 

s si ± = ~\M t \ 2 % st± ± IsflM^ejyj^i. (H3) 

This equation can be truncated by the method standardly used for moment expansion. One 
writes (v 2 i± ) — > u 2 si ± + a 2 i± and uses the equilibrium distributions for f S i± , s when evaluating 
the variance a 2 si ± = (v 2 i± ) — u 2 si± and the remaining integrals X psi ± appearing in the source 
term (113). 

Remarkably, expressions (112-113) show that the divergence of the axial current in fact 
corresponds to the first velocity moment of the kinetic equation for gQ d . Because a nontrivial 
CP- violation is tied to nontrivial complex phases in the pseudoscalar, or axial mass term, the 
fact that the source appears in the axial current nicely explains why the semiclassical source 
appears at first order in moment expansion in earlier semiclassical treatments [29, 15, 22]. 

The first source in (113) is to leading order in gradients spin-independent and does not 
violate CP. It is important for the phase transition dynamics however, in that it provides 
the dominant contribution to the friction on bubble walls from fermions [7, 9]. The second 
source is spin-dependent and CP-violating and it is thus responsible for baryogenesis. This 
is one of the main results of this paper, as it shows how the source from the semiclassical 
force enters to momentum integrated transport equations used in practical calculations. To 
promote equations (106) and (112) into transport equations for baryogenesis calculations 
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we still need to generalize them to include collisions as indicated in (46), which will be done 
elsewhere. 

5.1 Spontaneous source in relaxation time approximation 

We have so far shown how the source arising from the semiclassical force enters in the 
equations for currents, or equivalently momentum averaged transport equations. We shall 
now give a heuristic account on how sources have been modeled elsewhere in literature. 
The method very often used in EWBG considerations, apart from the WKB-computations, 
employs the relaxation time approximation for the kinetic equations. Here we discuss how 
the relaxation time approximation can be incorporated into our formalism, and then make a 
comparison with literature. Including collisions in the relaxation time approximation into the 
Liouville equation (72) results in the following kinetic equation for the distribution function 
f q for a charge q: 

(d t + v q ■ d s + F q ■ dj;)fq = - f JLll*L > (114) 

where r q = T^ 1 is the equilibration time for q (which we expect to be given by the rele- 
vant elastic scattering rate), F q the semiclassical force and f q o is the thermal equilibrium 
distribution function. In presence of a background field that violates q, one expects f q o to 
be shifted with respect to the naive thermal equilibrium, leading to a 'spontaneous' source 
that violates q. This source is more important for thick walls when the equilibrium f q pa f qQ 
is approximately attained on the wall. The spontaneous baryogenesis source was originally 
introduced by Cohen, Kaplan and Nelson [35, 16] in the context of two Higgs doublet mod- 
els, and then subsequently refined to include the m 2 -suppression in [41, 15]. The derivation 
was successively reconsidered in [36, 19, 20]. For example, in [20] the CP- violating vector 
current j q = f(dk z /(2n))fq for charginos in the MSSM was computed and inserted into the 
transport equations written in the relaxation-time approximation. 

The spontaneous baryogenesis source can be in our formalism obtained simply by inte- 
grating (114) over the momenta. The source then becomes the CP- violating contribution to 
the vector current (44), which to first order in gradients (or equivalently H) reads 

n Si - f si+ ~ j°_ - s\M t \^ / ^ 1 § (l + £(1 - f u )) , (115) 
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where we approximated f S i± by the equilibrium distribution function in plasma frame, f S i± — > 
l/( e ^/ T + l) « fu=F(s\M i \ 2 e' i /2u> 2 )df u ,/du>, where / w = l/(e w / T + l), and we used = 
ujduj. To make a comparison with literature, note first that the spontaneous source (115) is 
nonanalytic in |Mj| 2 . Since earlier attempts [36, 37, 17, 19] used expansions in powers of |Mj| 2 
to compute the spontaneous source, their results are at best incomplete. Consider next the 
CP-violating source for charginos in the MSSM. According to our equation (82) it is given 
by \M d \ 2 & = diag(m 2 h © / + ,m 2 _©'_), where m\Q' ± = ^(g 2 /A.)$f(fjt,m2)(h 1 h 2 )', showing the 
parametrical dependence (hih 2 )' on the higgs fields. This is in contrast with Refs. [19, 20], 
where a source proportional to hih' 2 — h 2 h[ was found and claimed to be important for 
baryogenesis. The origin of the difference may be in the fact that we made use of the 
constraint equations, which is necessary to obtain the correct results. 

Consider now the axial vector current. The corresponding spontaneous source can be 
easily obtained from (45): 

*jL± = E s kz j^f sd± = 0, (116) 

s k z — ± 

where we took f si ± — > l/(e w / T + 1). As a consequence, there is no spontaneous baryogenesis 
source from the axial vector current when computed in the relaxation time approximation. 

An attempt to compute the spontaneous baryogenesis source was made by Riotto [37], 
where the divergence of the vector current was computed in the Schwinger-Keldysh formal- 
ism [42] and then, based on [36], identified with a spontaneous source [17]. In this way he 
found an equation which formally reads: d^j^ ~ spontaneous source. According to Eq (106) 
however no source appears in the continuity equation for the vector current. Instead sources 
appear in the continuity equation for the axial vector current (110-112), which has not been 
so far considered in literature. 

6 Discussion and summary 

The question of a first principle derivation of CP-violating fluxes in transport equations has 
been the main theoretical challenge of recent work on electroweak baryogenesis. In this paper 
we derive the kinetic equations appropriate for EWBG in a systematic gradient expansion 
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starting from the (Dirac) equation of motion for the two-point Wightman function G < in the 
collisionless limit. The gradient expansion we use is well controlled and corresponds to an 
expansion in the de Broglie wave length divided by the wall width. In EWBG applications 
one typically has d-dB^w <C 1, so that such an expansion should be rapidly converging. 

We have shown that to first order in h the collisionless kinetic equations for both fermions 
and bosons can be recast as the Liouville equations for a single particle distribution function 
where the group velocity and the semiclassical force terms contain all quantum information, 
and in particular the CP-violating terms which source baryogenesis. These results agree with 
Ref . [22] , where the kinetic equations were obtained in the semiclassical WKB picture, orig- 
inally developed for EWBG problem in [29, 15]. The outstanding contribution of this paper 
is in a first principles derivation of these results in a completely controlled approximation 
scheme. We also derive the semiclassical force in the general case of N mixing fermions and 
in particular for the chargino sector in both the MSSM and NMSSM. Finally we prove that 
there is no CP-violating force at first order in h for scalar fields (102). Let us point out that 
the fact that the quasiparticle picture of plasma still holds to order Ti in gradient expansion is 
not surprising since the gradient correction for fermions from a pseudoscalar (CP-violating) 
mass condensate can be equivalently reformulated in terms of a 'classical' axial vector field 
condensate [29, 15]. 

We have also studied the vector and axial vector current equations, and showed that, 
while the vector current is conserved, the axial current contains both CP-conserving and 
CP- violating sources. This is to be expected, as CP-violation in particular is known to be 
caused by complex pseudoscalar (axial) mass terms. We have then pointed out that the 
axial current equation corresponds to the first velocity moment of the kinetic equation. This 
explains why the CP- violating source appears at first order in moment expansion [29, 15]. 
We finally made connection between the present results and literature where the continuity 
equations were written in the relaxation time approximation. In this context the source has 
been claimed to appear either as the vector current divergence [37, 17] or the time-component 
of the vector current j° [36, 19, 20]. However, we have shown here that the vector current 
continuity equation (106) in fact contains no source. We have also computed j°, and applied 
the result to the MSSM, and found a parametrically different result from Refs. [36, 19, 20]. 
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For simplicity here we consider only the collisionless limit in 1 + 1 dimensions. One can 
show [27] that generalization to the 3+1 dimensional case does not affect our discussions in 
any qualitative way. The question of how to consistently include collisions we postpone to a 
future publication. 
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